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We study the effect that initial state fluctuations have on final particle correlations in heavy ion collisions. More 
precisely, we focus on the propagation of initial perturbations on top of the expanding fireball using the con- 
formal solution derived by Gubser and Yarom for central collisions. For small perturbations, the hydrodynamic 
equations are solved by separation of variables and the solutions for different modes are added up to construct 
initial point-like perturbations, that are then allowed to evolve until freeze-out. The Cooper-Frye prescription 
is used to determine the final particle distribution. We present the two-particle correlation functions and their 
Fourier spectra obtained for different viscosities and different initial widths of the perturbation. We find that 
viscosity kills the higher harmonics, but that the Fourier spectra presents maxima and minima, similar to what 
is seen in the study of cosmic microwave background radiation. 



1. Introduction 

The main idea behind the present work was to study the evolution of an initial perturbation on top of the 
background produced in heavy ion collisions, and to calculate measurable quantities such as the two particle 
correlations and the Fourier coefficients of its expansion. However, before going into any of the details of our 
study we would like to start by talking about the background of the topic to be discussed and by giving some 
motivation for our calculations. 

When two relativistic heavy ions collide at experiments such as the Relativistic Heavy Ion Collider (RHIC) 
and the Large Hadron Collider (LHC) a new form of matter, the quark gluon plasma (QGP), is created. There 
is evidence [TH!] for the existence of a large elliptic flow -characterized by the Fourier coefficient v%- which is in 
agreement with hydrodynamical simulations, and thus supports the idea that the QGP behaves like a fluid and 
can be studied using hydrodynamics. For a long time hydrodynamic studies used smooth initial distributions 
from which the vi coefficient could be calculated but which due to symmetry considerations implied that all the 
odd coefficients in the Fourier expansion of the particle distributions were zero. The work by Alver and Roland 
[BJ, where initial fluctuations due to the random positions of the nucleons in the nuclei were included, showed 
that the V3 coefficient was of crucial importance for understanding the ridge and the double hump structure in 
two particle correlations. The v% and other large harmonics have now been measured experimentally [7] for 
a wide range of centralities and this together with the fact that the coefficient is present at even the most 
central collisions strengthens the importance of the understanding of the initial state and its fluctuations. 

In a previous work [5] we used the Glauber model to compute the initial deformations e n and angles ip n up 
to n — 6, and found that there existed correlations between all odd harmonics, n = 1,3,5. We understood 
these correlations to be evidence for the presence of hot or cold spots; initial perturbations to the background 
that will propagate during the dynamical evolution of the matter created in the collision. In this talk we will 
discuss the idea of sound propagating from these initial hot spots. The idea was introduced in [9] and followed 
up in 8 and ID] , and consists mainly in sound circles that propagate from the initial perturbation and expand 
until freeze-out reaching the "sound horizon" radius. The excess of matter due to these circles would give rise 
to a particle distribution with two peaks or horns which has also been observed in the calculations done in . 

This talk is mostly based on work done in |10j . where we studied the evolution of an initial perturbation on 
top of the flow developed in [TH H3] by S. Gubser and A. Yarom. To simplify our calculations and to gain 
understanding of the way in which perturbations propagate, and their effects on measurable quantities such as 
the two particle correlations, we study the case with only one Gaussian perturbation on top of the background 
temperature. Using the tools from [13j the initial state and its hydrodynamical evolution can be described 
analytically. It is important to notice that in order to do this we we will need to restrict ourselves to consider 
only central collisions with conformal matter satisfying e = 3p and thus with a fixed speed of sound during the 
whole evolution. We will also need to limit ourselves to small perturbations since we will be using linearised 
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solutions to the hydrodynamic equations. 

The plan that we will follow in this talk is to start by placing a Gaussian shaped initial perturbation on top 
of the background and allow it to evolve until freeze-out, where using the Cooper- Frye formula we obtain the 
particle distribution and from there the two particle correlation and the coefficients of its Fourier expansion can 
be calculated and then compared to experimental data. 



2. Perturbations on Top of the Fireball 

In this section we present the evolution of an initial perturbation on top of the background flow but, before 
going into this, it is necessary to understand the background, so we give a summary of Gubser's flow. 



2.1. Summary of Gubser's Flow 

Let us briefly discuss Gubser's flow which was developed in [T2] as a generalization to Bjorken flow [Uj and 
extended in |13j . The assumptions made by Bjorken are as follows: boost invariance along the beamline, rotation 
invariance in the azimuthal angle (j> around the beamline and translation invariance in the transverse plane. 
The first two of these premises are maintained by Gubser, but the third one is modified by a special conformal 
transformation in order to take into account the finite size and the transverse expansion of the fireball. The 
solution can be found for conformal matter with equation of state e = 3p ~ T 4 , which implies a speed of sound 
c s = In coordinates (r, r, (f>, r/) the four-velocity is given by 

Up = (— cosh k(t, r), sinh k(t, r), 0, 0) , (1) 

with 



/ 2(^ 2 T7* \ 

v± = tanh «(r,r) = 5 ^7 , (2) 

\ l + q 2 T 2 _|_ q 2 r 2 I 



while the energy density corresponds to 



«< 2 " 8,S (3) 



r±/3 (1 + 2q 2 (r 2 + r 2 ) + q±(r 2 - r 2 ) 2 ) 4/3 ' 



Here eo is just a normalization constant and q, which has units of inverse length, is the parameter that accounts 
for the finite size of the nucleus and for the radial flow. When q goes to zero the transverse velocity vanishes 
and Bjorken's solution is recovered. 

The solution that we have just described can be re-derived if one uses transformation from the r, r coordinates 
to a new set p, 9 given by: 

. 1 -q 2 r 2 + q 2 r 2 
smhp = — (4) 

2qr 

tan 6 = t^-t, 5-5-, (5) 

1 + q*r z — q z r z 

and re-scales the metric as follows 

ds 2 = r 2 ds 2 . (6) 

This was done in [T3] by Gubser and Yarom, and has the advantage that the fluid is now at rest (the only 
non-zero component of the four- velocity is u p = —1), so expressed in this new way the solution looks much 
simpler, and as we will see in the following section it is possible to analytically study perturbations on top of 
the background that it describes. The rescaled metric in the new coordinates reads: 

ds 2 = -dp 2 + cosh 2 p (de 2 + sin 2 6d<t> 2 ) + dry 2 , (7) 
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where, the new time coordinate is p and the new radial coordinate is 9. For the the case with shear viscosity rj 
the temperature in the rescaled frame is 

T ff sinh 3 p ^ fi 7 5 2 



T (cosh p)2/3 + 9(cosh p)2/3 2i?1 2 ' 6 ; 2 ' Smh P 

(8) 

where T is a normalization constant, H — rj/T 3 and 2-F1 is the hypergeometric function. 



2.2. Perturbations of the Flow 

In real heavy ion collisions the fireball is not smooth but it presents fluctuations due to the random positions 
of the nucleons in the nuclei. In order to take this into account, following the calculations done in |13| . we 
studied perturbations on top of Gubser's flow. Since we are interested in trying to understand the final two- 
particle correlations which present long range correlations in rapidity, we study only perturbations that are 
rapidity independent, this is perturbations in the transverse plane. The temperature and the four-velocity in 
the rescaled metric can be written as the background plus the perturbation in the form: 

f = f b (l + S) (9) 

^ = u hll + u lfl (10) 

where the perturbations correspond to 

S = 6(p,6,cj>) (11) 

ui„ = (O,w e (p,0,0), M0 (p,0,0),O) (12) 

By assuming that the perturbations are small, the hydrodynamic equations V^Tj,,, can be treated keeping only 
linear terms. By doing this one obtains a system of coupled first order differential equations, which in the 
ideal case may be decoupled to generate the following second order differential equation for the temperature 
perturbation 

d 2 S _L fd 2 S _L_cM 1 d 2 S \ 4 M 

^~3™^p\W + ^cW + ^W) + 3 t&nhp dp- - {U} 
Using separation of variables the solution is found to be 

5{p,e,4>) = 8i{p)Y lm {9,4>) (14) 

where the 8i(p) functions are given by combination of associated Legendre polynomials times cosh (p) 2 ^ 3 and 
the Yim are spherical harmonics. The perturbations to the velocity are as follows 

ue = vi(p)d e Yi m (9,<j)) (15) 
u<f, = Vi{p)d^Yi m {9,4>) (16) 



with the p-dependent part given by 



3cosh 2 pd^ 

Mp) = KiTiTdp- (17) 



In the viscous case the equations are more complicated, but separation of variables can still be used. The 9, <f> 
part of the solution does not change, but the p-dependent part does and it needs to be solved numerically 
in this case. The comparison between the ideal and the viscous case is presented in Fig{l] for three different 
values of 1. These plots show the oscillatory behaviour of the solutions, which is more evident for larger 
1; what begins as only a perturbation in temperature then becomes a perturbation in velocity and then 
keeps oscillating between them. The effect of viscosity can be seen for larger times, and it consists in 
a damping of the amplitude of the oscillation. This damping is greater for larger 1, so in the cases where 
viscosity is present the larger harmonics are suppressed, as we expected from the discussion in section II A of [S]. 
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Looking back to Eq |14[ we see that at any given time p the temperature perturbation is given in terms of 
spherical harmonics of the radial coordinate 9 and the azimuthal coordinate <fi, so we can use the completeness 
of these functions in order to write any initial condition as a sum of them by calculating the appropriate 
coefficients. In particular, at an initial time p = po we used a Gaussian perturbation in temperature with width 
s localized at 6 — 9q and <f> = <fio given by 



' (*-*o) 



cx e 



(18) 



and assumed zero initial perturbation to the velocity. By doing this we were able to write the temperature 
perturbation (in the rescaled frame) at any given time p as 



(19) 



where Ri{p) — Tf,( j o)(5;(p) and the c; m coefficients are calculated by using the orthogonality of the spherical 
harmonics and the initial conditions. In this way we may follow the evolution of the perturbation as it is shown 
in the three plots from Fig{2j What we find is that all the harmonics (we used up to I = 30) add up coherently 




r,o.s 



r,o.5t 
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Figure 2: Evolution of the temperature perturbation Ti(t, r, (j>) in the rescaled frame but in the regular coordinates for 
T = 1,4, 6 fm/c from left to right. (Taken from |10| ) 



to form the Gaussian initial perturbation, which then starts to propagate as an expanding circle on top of the 
background. 
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3. The Sound Circle and Final Particle Distributions 

In this section we present the single particle distribution, the two particle correlation and the power spectrum 
calculated from the perturbation that evolved on top of the background at freeze-out. 

3.1. Single Particle Distribution 

In order to calculate the particle distribution we used the Cooper-Frye formula [15] 

where the function / corresponds to the thermal distribution, S is the freeze-out surface and the overall minus 
is there because we work in the mostly plus metric. In this work we used an isothermal freeze-out, this is: the 
freeze-out surface was calculated from T{T,x,y,rj) = Tf which was solved to get E M = (rf(x,y),x,y,r)). The 
shape of the freeze-out surface, and thus the freeze-out time, depends then on the temperature distribution: 
where the temperature was higher, the freeze-out time will be larger because it will take more time to reach 
the temperature of the freeze-out. In this way the effects of the temperature perturbation are preserved, the 
expanding circle which had a higher temperature than the background will have a larger freeze-out time and it 
will be visible in the freeze-out surface. The intersection between the edge of the fireball and the circle give rise 
to the particle distribution. 
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(rad) <f> (rad) (rad) 

Figure 3: Single particle distribution (arbitrary units) for initial Gaussian perturbations with widths 0.4 fm (magenta) 
and 1 fm (blue), for three different viscosities n/s — 0, 0.08, 0.16 from left to right. 

The particle distributions as a function of the azimuthal angle <f> obtained for different initial widths of the 
perturbation and different viscosities are shown in Fig. [3] It is important to notice here that we have taken 
a qualitative approach; we are interested in the shapes of the distributions and in the effects of viscosity and 
initial width of the perturbation, so we used arbitrary units, which is why we have suppressed the labelling of 
the y axis. The main features of the distribution are two horns with a dip in the middle which are centred 
around the initial position of the perturbation. In the ideal case there appears a bump in the middle of the dip, 
but it becomes less important when viscosity is included. Other effects of viscosity are to damp and broaden 
the peaks and smooth out the curves. The width of the initial Gaussian perturbation also affects the final 
particle distribution; for a greater initial width the main characteristics of the curves are more pronounced: the 
maxima are larger and the minimum in the middle is more pronounced, also the width of the peaks is larger. 
This is because we kept the initial height of the perturbation and changed the width, so the total energy of the 
perturbation is larger, making the effects on the final distribution larger as well. 

3.2. Two Particle Correlations 

We would now like to qualitatively compare our results to some experimental data. In order to do so we 
calculate the two-particle correlation function by multiplying two single particle distributions dN /d<p from the 
previous section and averaging over the azimuthal position of the initial Gaussian perturbation. Since we are 
studying only perfectly central collisions this result depends only on the difference Aip — 4>i ~ 4>2- Our results 
are plotted together with data from ATLAS [TB] in Fig. [4] for comparison. We present two plots calculated for 
different viscosities and for an initial perturbation of width 1 fm (top of Fig. [4]) , where it is possible to see that 
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the main features of the two-particle correlation function resemble what is measured experimentally. There is 
one large central peak at A(f> — and centred around A</> = it there is a double hump, and since we did not 
include any rapidity dependence in our perturbation, these structures are elongated in rapidity. Although there 
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Figure 4: Top: Two particle correlations (arbitrary units) as a function of obtained for an initial Gaussian perturba- 
tion of width 0.7 fm for the ideal case (left) and for the case with n/s = 0.16 (right). Bottom: Two particle correlation 
(arbitrary units) as a function of A(f> obtained for an initial Gaussian perturbation of width 1 fm for the ideal case (left). 
Two particle correlation function measured by the ATLAS collaboration [16] for — 1% centrality (right). 

is a similarity between our results and the experiments, we must mention that they do not exactly match; in our 
case the flat zone centred at it is more extended than in the data. To understand this it is important to remem- 
ber that Gubser's flow that we use is an idealization, were matter is assumed to be conformal and to satisfy the 
equation of state e = 3p, which implies a constant speed of sound c s = 1/ \/3 throughout the whole evolution. 
This doesn't happen in real heavy ion collisions, where the speed of sound varies. The final radius of the sound 
circle will depend on the speed of sound and on the freeze-out time, thus the distance between the two bumps in 
the two particle correlation is also affected by these parameters. Since we use an idealized situation our results 
do not perfectly match the experiment, but it allows us to do a qualitative study and check that the picture of 
an expanding circle from an initial localized perturbation yields a similar result to the ones found in experiment. 

The top left plot in Fig. [4j for an initial Gaussian perturbation of width 0.7 fm and zero viscosity, 
presents some extra structure in comparison to the experimental data; next to the main peak at Ac/) = 
two smaller peaks appear. However, if viscosity is included (top right plot) or the width of the initial 
perturbation is increased (bottom left plot) these extra peaks disappear. As we will see below, both the initial 
width of the perturbation and the viscosity of matter affect large harmonics, which die for large values of 
these parameters, so we can conclude that the extra peaks are due to the large harmonics and that to match 
the experimental results the widths of the initial perturbations and/or the viscosity of the medium must be large. 
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The two particle correlation function can be expressed in terms of its Fourier expansion, given by 

dN 



dA<jy 



oc 



1 + 2 y% cos mA<p, 



(21) 



where the v„ are the flow coefficients in the Fourier expansion of the single particle distribution. Of course, if 
we now plot v„ vs n, we will be looking at the same information that is contained in the two particle correlation 
plots, but it is still interesting to look at the power spectrum and its features. In Fig.[5]it is possible to appreciate 
that the power spectrum appears to have a first maximum, then starts to decay and then secondary maxima 
appear. This resembles what is seen in the Cosmic Microwave Background (CMB) radiation spectrum where 
after a first maximum some secondary maxima follow (see, for example |17j). The plot on the left shows the 
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Figure 5: Spectral plots v\ vs n. Left: for viscosities n/s —0 (small dashed magenta), 0.08(dot-dashed red), 0.134 (solid 
green), 0.16 (dashed black) and width of the initial Gaussian perturbation of 0.4 fm. Center: for widths of the initial 
Gaussian perturbation of 0.4 (dot-dashed red), 1 (solid blue) fm and viscosity n/s = 0.08 . Right:from the ATLAS 
collaboration [15]. Both the left and center plots have been normalized such that the value of v^ for n = 3 matches the 
data. 



dependence of the Fourier coefficients on viscosity, and it can be seen that, as we mentioned earlier, viscosity kills 
the higher harmonics. The same effect is seen in the middle plot, this time we study the effect of the width of 
the initial Gaussian perturbation and we find that for larger initial widths the larger harmonics are suppressed 
in comparison to when smaller initial widths are used. This agrees with the study by Qin et al |18| . where 
it was shown how the larger widths of the initial Gaussian distributions suppressed the initial deformations 
e n with larger n. Our results are to be compared to the plot from ATLAS [19], which shows a large peak 
with a maximum at n = 3, the curve then falls, but at n = 9 a second maximum seems to appear. This 
second maximum is much smaller than than the first one, implying that a large width of the initial Gaussian 
perturbation and/or a large viscosity are necessary. 



4. Summary 

In this talk we presented the evolution of an initial state Gaussian deformation on top of the fireball, and 
the effects it has on measurable quantities such as the two particle correlations and the flow coefficients v n . 
Using the formalism developed in |13) by Gubser and Yarom, we were able to write the initial perturbation as 
a sum of the solutions to the hydrodynamic equations, and to see the evolution of the perturbation in time as 
it propagated from its center as an expanding circle, that by the freeze-out time had reached the edge of the 
expanding matter, and contributed two horns to the single particle distribution. 

We looked at the two particle correlations and to the flow coefficients calculated for different widths of the 
initial Gaussian perturbation and for different viscosities. The variation of these two parameters proved to 
have similar results, the larger their values the smoother the particle distributions, and the more suppressed 
the values of the higher harmonics become. Comparing our results for the spectral plots with the data we find 
that for them to agree a large initial width of the perturbation and/or al large value of viscosity would be 
necessary, because the magnitude of the second peak is much smaller that that of the first. 
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It is necessary to mention that we worked with an idealized situation, where we looked only at the case 
with only one initial Gaussian deformation but in reality the initial state fluctuations are determined by the 
random positions of the nucleons in the nuclei and many such deformations are expected. In order to take this 
into account it would be necessary to average over many initial deformations using probability distributions to 
calculate their locations and amplitudes. Other assumptions that we used that do not hold in real heavy ion 
collisions are that we took the matter to be conformal throughout the process, and the initial perturbations 
to be small, in order to be able to use the analytical tools at hand. Still we find it remarkable, that even 
though we treated a very idealized situation in which we assumed that all harmonics add up coherently to form 
the expanding circle from the perturbation, the results that we obtained agree qualitatively quite well with 
experimental data. 
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